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Abstract. Let o-{t,t') be the sigma-algebra generated by the differences Xs —X^' 
with s,s' £ where {Xt)~ao<t<oo is the fractional Brownian motion with 

- - - Hurst index H £ (0, 1). We prove that for any two distinct timepoints ti and t2 

, the sigma-algebras a(t-^_^^ti+s:) and aft2-s,t2+£) are asymptotically independent 

' as £ \ 0. We show this in the strong sense that Shannon's mutual information 

^SJ , between the two cr-algebras tends to zero as £ \ 0. Some generalizations and 

^ ' quantitative estimates are provided also. 

o 

:z; 

^ ' 1. Introduction 

Let X = (Xt) — oo<t<oo be the standard fractional Brownian motion (FBM) with 
P_( ■ Hurst index H G (0,1). Thus, Xt is a centered Gaussian process with stationary 

(~| , increments and variance function EXf = \t\'^^ (see, e.g., [7, 11]). The parameter 

I value H = yields the standard Brownian motion. FBM is a ii-self-similar pro- 

cess, that is, {Xat)^={a^ Xt), where = stands for the equality of finite-dimensional 
distributions. 

When H ^, the increments of X on disjoint time intervals are always dependent 
Q>^ ■ — negatively correlated for i/ < ^ and positively correlated for if > |. Moreover, 

O ■ when H > ^, the sequence Xi, X2 — Xi, X^ — X2, ... is long-range dependent, i.e. 

Y,°Z^B Xi{Xi+i - Xi) = 00 (see [2, 11]). FBMs with if > i are often used in 
applications as a mathematical model for far-reaching dependence. 

However, as we show in this paper, 'small and distant' events in FBMs are nev- 
ertheless asymptotically independent. This holds both as asymptotic orthogonality 
of the Gaussian subspaces generated by the processes {Xt)\t\<i and (Xn+t — -'^n)|t|<i 
as n 00, and in the stronger sense that the mutual (Shannon) information 
^ I I{{Xt)\t\<:i '■ {Xn+t — -'^ra)|t|<i) is finite and decays to zero as n ^ 00. By self- 

' similarity, this is equivalent to considering the increment processes around two fixed 

timepoints, {Xg+u — Xs)\u\<e and {Xt+u — Xt)\u\<ej as e \ 0. We propose to call 
this latter property local independence. 

Our paper was motivated by [8], where FBM's local independence property was 
needed, but attempts to find this result from literature were unsuccessful. Very 
recently, however, J. Picard [9] has proven the asymptotic orthogonality result using 
a different technique. The more functional analytic approach of the present note 
has the advantage of giving very precise estimates both for the rate of asymptotic 
orthogonality, and for the much stronger property of asymptotically vanishing mutual 
information. 
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The structure of the paper is as follows: in the first section we briefly recall certain 
facts about Sobolev spaces with fractional smoothness index - these spaces are the 
main tool in our approach. We have tried to make the exposition readable for the 
readers with no previous knowledge on these spaces. The second section reviews the 
basic facts on the Gelfand-Yaglom theory of mutual information between Gaussian 
spaces. The third section contains the proof of our main results. The results are 
obtained in a quantitative form in terms of the relative size of the time intervals 
involved. Finally, the fourth section briefly considers the higher dimensional case 
and states open questions. 

2. Preliminaries I: the fractional Sobolev spaces 

We shall apply the common notation for uninteresting constants. They will all be 
denoted by the letter c, and its value can vary inside a single estimate. The notation 
a ^ b means that the ratio of the (positive) quantities a and b stays bounded from 
below and above as the parameters of interest vary. The inner product of elements 
(f) and f/j of a Hilbert space H will be denoted as ((p,xp)-n, and the angle <{A,B) 
between subspaces A and B of 7^ is defined by 

cos«(^, B)) := sup { -.UeAVeB 

Suitable references for this section are e.g. [10, Section 6] or selected parts of [12]. 
Vastly more information can be found in Triebel's monographs, like [15]. Actually 
only very little of the theory of Sobolev spaces is needed, and we try to be as self- 
contained as possible. 

The Fourier transform of a tempered distribution / on R" is defined as 

/(O := (27r)-"/2 / e-'-<f{x)dx. 

We shall employ the notation (A, n) for the distributional pairing, assuming that it 
is well-defined for A and /i. Recall that the convolution A * is always defined if A is 
a Schwartz distribution and (p G C^(R"), and its Fourier transform is the product 
(27r)"/^^ A. Moreover, by the definition of the Fourier transform, the Parseval identity 
can be written in the form _ 

(A,0) = (A,0). 

Let s G R. The Sobolev space Ty^'2(R") is defined as the Hilbert space of tem- 
pered distributions / on R" such that the Fourier transform /(^) is a locally inte- 
grable function with the property 

(1) ||/||.,2 := WfWws,^ := (/^^ + l^l'r) < oo. 

Our normalization constant for the Fourier transform makes sure that VF*''^(R") = 
L^(R"') isometrically. 

In the distributional pairing, the isometric dual of W^'^iR") is W-'''^{W). More- 
over, the norm increases as s increases, and for integers € N we have that 

(2) 11/111,2 ~ / {\f{x)\'+ J2 \f^''H^)f)dx. 

\a\=k 

Obviously all these spaces are translation invariant, and one may verify that multi- 
plication by an element in Co"(R'*) is continuous. 
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We next recall the homogeneous Sobolev spaces M^*''^(R"').The norm is replaced 

by 

(3) ll/llvK..(R„) := (^^ mWd^y < oo. 

This norm is certainly well-defined at least for all / G C^(R"), although even then 
it may take the value oo if s < — n/2. In defining the Hilbert space W^*'^(R") there 
indeed arises some complications in the definition, due to the fact that the l^p* can 
be either 'too big' or 'too small' near origin. However, for our main result it is enough 
to consider the case n = 1 and \s\ < 1/2, and then these difficulties disappear. For 
these values of the parameters the homogeneous spaces are simpy defined as the 
(inverse) Fourier transform of the weighted space L^(R), where the weight is of the 
form n{d£,) = jCP'^- I^y Cauchy-Schwartz any function in this weighted space is a 
locally integrable function, and thus defines a distribution in a natural way. On the 
other hand, every Schwartz test function belongs to this weighted space, which can 
be used to show that Cq°(R) C VF*'^ is a dense subset. Moreover, the isometric 
duahty (^"^'^(R))* = W'^''^{'R) holds via the distributional duality 

= / 4'(x)ip{x) dx. 

The pairing is originally defined only for test functions, but it extends to elements 
e VF'^'^(R) and ijj G VF"*'^(R) by continuity and density. 

We then fix s G (—1/2, 1/2) together with an open interval / C R (/ can well be 
unbounded) and define the Sobolev-functions over this interval. First of all we denote 
by Wq'\I) the closure of C^{I) in the space W*'^(R). Clearly all the elements in 
Wq' (/) are distributions supported on /. We will also need the space W'^''^{I) which 
consists of restrictions of elements of W^''^{'R) on the interval /. Thus VF*''^(/) = 
{9\i '■ 9 ^ VF*'^(R)}. This space is naturally normed by the induced quotient norm 

■= iiif{ll5llw?.,2(R) : 9\i = /}• 

In a similar vain one defines the non-homogeneous space W^''^{I) by setting W^''^{I) = 
{9\l '■ 9 € W^*'^(R)} and introducing the quotient norm 

:= inf{||£'llw''.2(R) : 9\l = /}• 

This definition makes sense for all s G R. One may easily verify that ||/||^i,2(/) ~ 
Iiif^i^) + f'^{^))dx, where /' is the distributional derivative of /. 

bmce Wo \l) C l'F-^'2(R) is a (closed) subspace, we deduce by standard Hilbert 
space theory that isometrically 

(4) {W^'^{I)y = W-''^{I) and {W'''^{I)y = W^"'^{I) 

through the pairing {(f), ip) = Jj (p{x)^{x) dx (extended again by continuity). There 
is thus a natural isometry G : W~^'^{I) — ^ Wq''^{I) in such a way that 

(5) ((f), Gip)^.,2^j^ = (f){x)ij{x) dx 

for smooth elements <p and ijj. Again this extends for any G Wq'^{I) and if) G 
Ty~*'^(/) by continuity. 

In the Lemma below the assumption \s\ < | is crucial. 
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Lemma 1. Let s G (~ 2' 2) let I C H be an open interval of length 1. 

(i) Multiplication by the signum function extends to a bounded linear operator on 
VF*'2(R). In other words, ||x(-oo,o)/ll^ys,2(R), ||X(o,oo)/IIm7«,2(r) < c\\f\\^,^2^^^) for all 
f G Co°(R). The same statement remains true i/T^*'^(R) is replaced by Ty*'^(R). 

(ii) W^'\l) = {fe W^'^iK) : supp(/) C 7}. 

(iii) We have Wq'\I) = W''^{I) = W'^il) with equivalent norms (the constant of 
isomorphism does not depend on the location of the interval I). 

(iv) There is a continuous inclusion W^''^{I) C Wq''^{I), and this natural imbedding 
is a Hilbert- Schmidt operator. 

Proof, (i) The statement is well-known, see [15, First Lemma in Section 2.10.2.]. 
Actually, up to a constant the multiplication by the signum function corresponds 
to the action of the Hilbert transfrom on the Fourier side. Hence the claim follows 
from the fact that is a Muckenhoupt ^^-weight on R for any s G (—1/2, 1/2), 
see [13, Corollary, V.4.2, V.6.6.4]. In a similar way, by checking that (1 + \£,\'^Y is a 
Muckenhoupt weight one obtains the statement concerning W^'^. 

(ii) Let / G VF'''^(R) with supp(/) C /. We will show that one may approximate 
/ in norm by the elements of C^{I). The dilation A — > /(A-) is a continuous map 
from a neighbourhood of 1 into 1^*'^(R). Hence, by approximating / with a suitable 
dilation we may assume that supp(/) is contained in I. Finally, we then obtain the 
required approximant by a standard mollification. 

(iii) By the translation invariance of the spaces, the independence on the location 
of the interval / is obvioTis. The first equality is an easy consequence of parts (i) and 
(ii). Towards the second equality, let us first verify that C^{I) is dense in VF^'^(I). 
By part (i), if / G VF*'^(/) then also xif £ P^*'^(R), where xif stands for the zero 
continuation of / to R. Exactly as in part (ii) we show by dilation and convolution 
approximation that xif is in the closure of Cq°(/) in PF*'^(R), which clearly yields 
the claim. 

Hence it remains to show that 

(6) \\f\\w^:^(n) - \\f\\w^,^n) / G C^{I). 

We may clearly assume that / = (0, 1). Let us first consider the inequality 

(7) ll/llw».2(R) < c||/||^^,,.(j^). 

This is immediate if s < 0. If s G (0, 1/2) we choose a cut-off function (p G Cq°{—1, 2) 
such that = 1 on the interval [—1/2,3/2]. Let us decompose 

where fi = X[-i,i]f, and ^ = /-/i- Then obviously ||0/2||i^^,2(r) < c||/2||iy.,2(R) < 
c|l/llw^.,2(R)- Moreover, 

/i(x) = ^|_\-«/(Ode, f[{x) = -^f^e^<m)d^ , 

where, by Cauchy-Schwarz, j\ \f{0\d^ < c||/||p^?,.2(R)- Hence ||/i||oo + ||/(||oo < 
c||/||^y.,2(R) and we obtain that ||(/'/i||Ty«,2(R) < ||<?!>/i||w^i,2(r) < c||/i||^;j^,,2(r)- By 
combining these estimates (7) follows. 
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In turn, the converse inequality 
(8) \\f\\w^,^R) < 4f\\w^,Hn)- 

is immediate if s > 0. It clearly follows for negative s G (—1/2,0) if we verify that 
in our situation ||/||l<x>(_i^i) < c||/||^ys,2(j{^). This is seen by observing that 

/(e) = -^(/(x),</,(x)e-^-), 



where sup_i<^<i \\(t>{x)e '^^||w^-.,2(r) < c. 

(iv) By part (iii), the claim is a consequence of the well-known Hilbert-Schmidt 
property of the inclusion W^'^{I) C W^''^{I). Since we have not been able to find 
a convenient reference, the simple proof is sketched here. We may assume that 
/ = ( — 1/2,1/2) so that / C (— 7r,7r] =: T, where T stands for the 1-dimensional 
torus. By applying a simple extension one may consider the spaces in question as 
closed subspaces of the corresponding Sobolev spaces H^{T) and H^(T) on the torus, 
where for / = ^^=-00 «ne'"^ and u e R one sets = Er=-oo(l + l"l)'"l«nP 

(see e.g. [10]). By considering the natural orthogonal basis ((1 + |n|)~"e'"-^)5^_o^ we 
see that the embedding H^{T) C H^(T) is equivalent to the diagonal operator with 
the diagonal elements ((1 + \n\y~^)^^_^. This is Hilbert-Schmidt as Yl'^=-ooi^ + 

□ 

We shall need the formula for the Fourier transform of the function Uaix) := 
where x G and a G {0,n). It is well-known, see e.g. [12, V 1. Lemma 2, p. 117], 
that 

(9) n«(6 = dn,a|C|"-", where := 2"/2-"Ii^^-^. 
Lemma 2. Assume that s G (—1/2, 1/2). 

(i) Let a > 0, a 1 and denote fa{x) = (1 + |a;|)"". Then fa G W-%K) for 
a > 1/2- s. 

(ii) Let a > 1/2 + s, a ^ 1. Then for any k > there is a constant c{a,s) > 
such that \\{k — ■)~"llvi?-s,2((_oo 0)) ~ c{a, s)k^^^'^^~'^ . In other words, 

(10) sup / {k-x)-''(p{x)dx = c(a,s)fe^/2+3-a_ 

Proof, (i) Choose a smooth cut-off function (p G Co°(R) such that ^ = 1 in a 
neighbourhood of the origin. Compose 

fa{x) = <^(x)A(x)+(l-(/)(a;))(/a(x)-|x|-")+(l-<^(a;))|x|-" =: gi{x)+g2ix)+g3{x). 

Obviously 51 G L^{R) n VF^'2(R) C W^''^(R) for all |s| < 1/2. An easy eastimate 
shows that the same holds for g2- Moreover, we observe that {d/dx)g'^ G L^(R). 
Hence |CPl53(0P < 1- Thus the inclusion 53 G 1^*'^(R) holds if and only if the 
integral/^, |?Pl?3(0P is finite. 

Consider first the case a > 1. Then 33 G L^(R), so that 53 is bounded and 53 G 
W*'2(R) for all \s\ < 1/2. Assume then that a G (0, 1). Then gsix) - \x\-°' G -L^(R), 
so that (9) yields [53(6 - di,a\^\'^-'\ < C. Thus {^^3(0? < 00 exactly for 
s > 1/2 -a. 
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(ii) The definition of the homogeneous Sobolev norm yields the scaling rule 
ll<?^(^")llvKs,2(R) = k^~^^'^\\(j)i-)\\y^^s,-2(jiy By using this fact, duality, and a substitu- 
tion X = ky in the integral we are reduced to showing that 



r {i+\x\ 

J —oo 



sup / (1 + (p{x)dx <oo. 
0eCo°°(-oo,o) 

By duality this follows immediately from the fact that (1 + G W~^''^{'R) 

according to part (i) of the Lemma. □ 

We finally remark that all the results stated in this section remain valid with 
identical proofs for the Sobolev spaces that contain only real- valued functions. 

3. Preliminaries II: mutual information between Gaussian subspaces 

In this section we present the needed facts from the Gelfand-Yaglom theory of 
mutual information between Gaussian subspaces. In order to recall the general con- 
cept of mutual information, let {ft, T , P) be a probability space, and let A and B be 
sub-cj-algebras of JT. The mutual (Shannon) information between A and B is defined 
as [4] 

Here the supremum is taken over all v4-measurable partitions O = Uj=i ^fc ^^'^ ^~ 
measurable partitions Q = IJ^i -Bfe of the probability space (n,m > 1, P{Aj) > 
and P(Bfc) > for all j,k). 

For random variables X : Q ^ E and Y : Q ^ F, where E, F are measurable 
spaces, we set I{X : Y) := I{a{X) : ij{Y)). Let (resp. fiy-, l-J'(x.Y)) be the 
distribution (measure) of X (resp. Y, {X,Y)) in the space E (resp. F, E x F). 
Then, one may check that I{X : y) = oo if the measure IJ^[x,Y) absolutely 
continuous with respect to the product measure fix ® fJ-v- Moreover, in the case 
where tJ-(x,Y) denote p = and have the formula 



(11) I{X:Y)= [ log{p)d{fix^l^Y). 

JxxY 



The KuUback-Leibler information characterizes the shift from a probabiUty mea- 
sure /X to another probability measure u on the same measurable space, and it is 
defined as 

oo, otherwise. 
Shannon's mutual information can be expressed in terms of the Kullback-Leibler 
information as 



(12) I{A : B) = lKL{PiA,B) ■■ Pa ® Pb), 

where P(a,b) denotes the unique probability measure on (17 x ^ x B) satisfying 
P{A,B}{^ X B) = P{A n B) for A e A, B e B. Actually, this is obtained from (11) 
by letting X (resp. Y) be the identity map {fl,A) (resp. the identity map 

in,j^)^{n,B)). 

The following properties of mutual information are most conveniently proven 
through the relation (12). 
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Theorem 1. (i) I[A ■ B) > and equality holds if and only if A and B are inde- 
pendent. 

(ii) I{A : B) is nan- decreasing with respect to A and B. 

(iii) If An ] A and Bn T B, then I{An : &„) T I{A : B). 

(iv) // An i A and Bn [ B, and if I{An ■ Bn) < oo for some n, then I{An '■ Bn) i 
I{A:B). 

When X and Y are finite- dimensional random vectors such that {X, Y) is a non- 
degenerate and centered multivariate Gaussian, one may easily compute by using 
(11) that 



2 " det(r(x,y)) 

where Tz denotes the covariance matrix of a Gaussian vector Z. In particular, the 
information between random variables X and Y with bivariate centered Gaussian 
distribution is 

(13) I{(^{X) : (t{Y)) = -logsin<(X,y). 

The theory of Shannon information between Gaussian processes was developed by 
Gel'fand and Yaglom [5]. Their fundamental discovery was that one may express the 
information between two closed subspaces A and S of a Gaussian space Q in terms 
of the spectral properties of the operator T := PaPbPa, where Pa and Pb stand for 
the orthogonal projections on A and respectively. In order to explain their result, 
and for later purposes, we first recall some basic notions of operator theory. 

Let S : £^ — > F be a bounded linear operator between the separable Hilbert spaces 
E and F. Let {ei}ig/ be an orthonormal basis for E. The Hilbert-Schmidt norm of 
S is defined as 

\\S\\hs{e,f) '■= ll'S'eillf)^^^. 

This definition does not depend on the particular orthonormal basis used. In case 
< oo we say that is a Hilbert-Schmidt operator. Also it is clear that if E 
(resp. F) is a Hilbert subspace of a larger space E (resp. F), then ||S'Pg||^^^^ = 
l|5'llfl5(£;,F)- In this sense it is not important to keep exact track on the domain of 
definition and image spaces, and one usually abbreviates ||'S'||h5(e,f) = ||'S'||_ff5. For 
products of bounded linear operators between (perhaps different) Hilbert spaces we 
have 

(14) \\TS\\hs<\\T\\hs\\S\\ and ||.Sr||^^5 < ||5||||r||^^5. 

Let us then assume, in addition, that S : E ^ E is self-adjoint and positive 
semi-definite, S* = S and S > 0. Then one may always define the trace of S by 
setting 

tv{S) ■.= J2iei,Sei) 

iei 

Thus, tr (S") G [0,oo]. In the case that tr (5) < oo we say that S is of trace class. 
Every trace class operator S is compact, and since we also assume > 0, it has a 
decreasing sequence of positive eigenvalues Ai > A2 > • • • > 0, where each eigenvalue 
is counted according to its multiplicity. It follows that 

(15) tr(5)=^Afe. 

Afc>0 
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We finally observe that 'd S : E ^ F is any bounded linear operator, tlien S*S > 
is self-adjoint, and we may compute 

(16) tv{S*S) = \\S\\ls- 

Let us then go back to the situation where A, B are closed subspaces of a Gaussian 
Hilbert space Q and state the result of Gelfand and Yaglom. Again Pa and Pb 
stand for the orthogonal projections to the subspaces A and B, respectively, and 
I{A : B) := I{a{X : X e A} : a{Y : Y G B}). 

Theorem 2. [5] Denote T := PaPbPa- The mutual information I{A : B) is finite 
if and only if \\T\\ < 1 {i.e. <{A,B) > 0) and the operator T is of trace class. 
Moreover, in this case 

(17) i{A:B) = ^ Yl ^^g^rrr^' 

k■.\^.>0 ^ 

where Ai > A2 > . . . are the eigenvalues ofT in the decreasing order repeated accord- 
ing to their multiplicities. 

A nice sketch of the derivation of the formula (17) is included in a form of ex- 
ercises in [3, pp. 68-69]. Assume that T is of trace class, and let ^1,^2,... be 
an orthonormal basis of TQ consisting of eigenvectors corresponding to the non-zero 
eigenvalues Ai > A2 > . . .. It is not difficult to see that {PsZi} is an orthogonal basis 
of PbPaPbQ , and, moreover, these bases are mutually orthogonal: {Zi,PBZj) = 
for i ^ j. Since orthogonality implies independence in tlie case of Gaussian random 
variables, it follows that the information between (t{A) and cr{B) can be expressed 
as the sum of the informations within the pairs {Zi, PbZ^), given in (13): 

I{A:B) = -^logsin<(Z,,PB^i) 

i 

i i 

Note that since <{A,B) = mii <i{Zi., PbZi) , the information between subspaces can 
be infinite even when they have a positive angle. 

By invoking the Taylor series of x 1— > log(l/(l — x)) we obtain for x G [0, 1) that 



(18) X < log(— ) = E £ - + 2-^(13^) S + 

k=l ^ ' 

Observe also that T = {PbPaT{PbPa) and ||r|| = WPbPaW^- Moreover, Ai < 
II-Pb-PaII < II^b-PaIIhs- By combining these observations and the facts (15)-(18) we 
obtain a formulation suitable for our purposes: 

Corollary 1. The angle between the spaces A and B satisfies cos{<{A, B)) = \\PbPa\\ 
We have I{A : B) < 00 if and only if \\PbPa\\ < 1 and \\PbPa\\hs < 00. Moreover, 
in this case 

(19) \\\PbPa\\Is < HA-.B) < + 

Observe that the above estimate is asymptotically precise in the limit ||PbP^|| 
0, or, equivalently, as <{A, B) — > 7r/2. Especially this is true in the limit I{A : B) 
0. 
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4. Statement and proof of the main results 

In this section we consider the asymptotic independence of the local spaces of 
FBMs. To be more exact, let us first define for any set 5 C R 

Es := spEn.{Xu — Xy : u,v G S}, 

and the shorthand notation 

We consider the following two notions of local independence. 

Definition 1. We say that the stochastic process X possesses local independence in 
the weak sense, if for any distinct ti, t2 

We say that the stochastic process X possesses local independence (in the strong 
sense), if for any distinct ti, t2 

The term 'weak' corresponds to its use in 'stationarity in the weak sense'. 

We will consider integrals of the form j-^Xt(t){t)dt for smooth and compactly 
supported functions (p. The definition of the integral poses no problems since t^Xf 
is continuous with respect to L^-norm of random variables, whence it can be e.g. 
defined as the limit of the corresponding Riemann sums (or as a Bochner integral). 
Let us start with two simple lemmata. 

Lemma 3. For any T G R and a > the elements 

(20) [ (t>'{t)Xtdt, (f)eC^{T,T + a) 

are dense in E(^x,T+a)- 

Proof. By observing that j-^(t)'{i)Xtdt = j-^cp' {t){Xt — Xj^^,^/2) d,t we see that the 
elements in question are contained in ET,a- Conversely, let (p ^ C'o°(R) satisfy 
Jj^(/)(t) = 1. Denote ^^(x) = e^^(j){xe). By the L^-continuity we have that for any 

ti,t2 e (T,r + a) 

Xt, - Xt^ = lim( / Xu{(i)e{ti +u)- (t)e{t2 + u)) du. 

Observe that we may write ■0' = ^e{ti + •) — 4>e{i'2 + •) for suitable ip G Co°(R). This 
yields the claim. □ 

Next we verify that the L^-norm of a random variable of the form (20) equals the 
norm of in a corresponding homogeneous Sobolev space. For later purposes we 
first state an auxiliary result that is valid in all dimensions. 

Lemma 4. Assume that H G (0,1) and the functions E Cq°{W') satisfy 

/j^n 4>dx = ijj dx = 0. Then 

(21) / / l{\u\'^^ + \v\'^^ -\u-v\'^")^{u)^{v)dudv 
7r" 7r" 2 

= -2 TT / T{-H) K'?.V')ty-n/2-if.2(R„)- 
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Proof. We first claim tliat for a £ (0, n) 

(22) / / \u - v\-''4>iu)i^ dudv = i27rr dn,a [ \^r-''mW)d^- 

JR" JR" JR" 

This is immediate by (9) and the Parseval formula since the left hand side above 
can be written as g'tp dx where g is obtained as the convolution g = Ua * (p, 
whence its Fourier transform equals (27r)"/^d„.Q,|,^|"'~"0(,^). By the assumption we 
see that the Fourier transforms of (p and V' satisfy |^(!>(0|, [tpiOl ^ near the origin. 
Moreover, they decay polynomially as |^| — > oo. These observations verify that the 
right hand side of (22) is analytic as a function of a in a neighbourhood of the open 
line segment a G (— 2,n). Since the left hand side of (22) is likewise analytic in the 
same neighbourhood we deduce by analytic continuation that (22) holds true for all 
a G (— 2,n). The claim follows as we substitute a = —2H in (22) and observe that 
by Pubini the terms \u\'^^ and \v\^^ make no contribution to the integral in the left 
hand side of (21). □ 

Corollary 2. Let H G (0, 1) and assume that ^i, (/'2 G Cq°(R) are real-valued. Then 
E (^{J^4>[it)Xtdt){J^4>2{t)Xtdt)^ = aH{cl>i,h)^i_^,., 

where an '■= sin(7ri?)r(l + 2H) > 0. E,specially, there is an isometric and bijective 
isomorphism J : -E(_oo,oo) ~^ VF^/^~^'^(R) so that for each interval {t,t') C R we 
have JiE^t,t')) = W^^''~'''\it,t')). 

Proof. Let us denote 



^ := E (^i j^ ^[it)Xt dt){ j^ (t)'2{t)Xt dt) 



By the definition of the fractional Brownian motion with the Hurst parameter H G 

(0, 1) we have 

(23) A = ]- f (l)'^{u)24>'{s){\s\^^ + \u\^^ -\s-u\^^)dsdu 

2 JrxR 

(24) = ai7(0i,02)H?-l/2-H,2 = dif (01, 02)^1 -H,2- 

Above we used Lemma 4 to obtain the first equality. Observe that the functions ^'^ 
and (p^ automatically have mean zero. The last equality follows directly from the fact 
that the Fourier transfrom of <^'j equals i4(j)j{S,), j = 1,2. The constant is simplified 
by applying the standard formulas for the Gamma functions, see e.g. [1, 5.2.4]. The 
last statement of the Corollary follows immediately by Lemma 3. □ 

Remark 3. Note that an takes the value 1 for iJ = 1/2 and tends to zero as — > 1~ 
or H 0+. 

Let us observe that if the supports of (pi and 4>2 are disjoint, we are free to integrate 
by parts in (23) and obtain the formula 

(25) E ("(^ ^\ (t)X, dt){J^ 4>'2(t)X, dt)) 

JRxR F ^1 

Here it is interesting to observe the sign of the factor H(2H — 1) for different values 
of the Hurst parameter H. 
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We are now ready to prove the main result of the paper. 

Theorem 4. Fractional Brownian motions with H G (0, 1) possess local indepen- 
dence. Moreover, there is a constant rj^ > (with r// > for H ^ 1/2) such 
that 

cos{<{Et,,e, Et,,s) = rnie/lti - tsl)'"'^ + 0{£^-^") as e ^ 0, 
and (with some 6h > 0) 

Proof. By scaling invariance and stationarity it is equivalent to show that 

(26) cos(<(£;(o,i),£^(fc,fc+i))) = rnk'^"''^ + 0{k^"'^) as k ^ oo and 

(27) /(i?(o,i) : i?(fc,fc+i)) = ir|,fe4^-^ + 0(^4^-5) as A; ^ oo. 



Denote s := 1/2 - H e (-1/2,1/2) together with A := H^o'^(/c,A; + 1) and B := 
Wq'^{0, 1), considered as subspaces of the Hilbert space W''^'^{'R.). Let Pa (resp. Pb) 
stand for the orthogonal projection on A (resp. B). We will consider the operator 

S := Pb ■■ A^ B. 

Since S = {PbPa)\a and {PbPa)\a^ — 0) obtain that = ||PbPa|| and 
||5'||_f/5 = ||-Pb-Pa||h5- Hence Corollaries 1 and 2 yield that 

(28) cos(<(£;(o,i),^(ik,fe+i))) = \\S\\ 
and 

(29) \\\S\\ls < /(%i),^(fc,fe+i)) < ^11^111^5(1 + ll^ll) 

as soon as \\S\\ < 1/2. 

In order to estimate the norm and the Hilbert-Schmidt norm of the operator S 
we will make use of the decay of the kernel in (25), and the even faster decay of its 
derivatives. For that end we need to first factorize S properly through a suitable 
integral operator. Assume thus that k > 2 and (f) G CQ°{k,k + 1) G A. Then by 
definitions and formula (25) we see that Scf) E B is the unique element that satisfies 
for each V e C^(0,1) 

f d)(v]ih(x] 
(Sep, ^)^.2(o,,) = (</>, i^)^s..^^^ = H{2H - 1) j^^^^^^^^^^^^^ 1^ _ y\2-2H d^dy 

(30) = f ^{x){R(l)){x)dx, 
where R stands for the integral operator 

R(t){x) := H{2H - 1) f ^^^^ dy. 

J{k,k+i) \x - yV'^" 

By the smoothness of the kernel we immediately see that R is well-defined and, in 
fact 

^ Ji {Wf(k, k + l))G W^'\{<d, 1)). 

Let G : W~^'^{0, 1) — > VFq'^((0, 1)) be the isometric isomorphism from (5). According 
to (30) we may factorize 

S = GR. 
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Let V : Wo'^(A:, k + l)^ VF^^'2(0, 1) stand for the one-dimensional operator 
V(t>{x) := I <P{y)dy, for (0,1). 

•J(k,k+1) 

Thus V(j) is constant on (0, 1). We decompose 

S = H{2H -l)k'^"-^GV + G{R-H{2H -l)k^"-^V). 
If we show that 

(31) \\(^R-H{2H-l)k^^-^v) ■.wf{k,k + l)^W-'^\k,k + l)\\HS 
= Oik^''-^), 

then, according to (28)- (29) and the fact that for the one-dimensional operator GV it 
holds that IIGV^IIijs = (the value is independent of k), both of the asymptotics 

in (26) follow immediately. Here we also keep in mind that the Hilbert-Schmidt norm 
always dominates the operator norm. 

Observe towards (31) that for x G (0, 1) and (f) G GQ^{k, A; -|- 1) we may write 

{{R - {H{2H - I)fc2^-V)</.) {x)=c f u{x, y)<j){y) dy, 

J{k,k+l) 

where a simple computation shows that the the kernel u{x, y) = \x — y\'^^~^ — k^^~^ 
satisfies 

||((^)"(^)M(^'-)IU-(M+i)<cfc2^"', a,/3e{0,l}, XG(0,1). 

By Lemma l(iii) we have || • \\w-'>,2(^(^k,k+i)) - H ' \\w^-^{ik,k+i))- Hence the previous 
estimates yield for fixed x G (0, 1) the estimate 

(32) \\u{x, •)llvK-.,2((fc,fc+i)) < Mx, Oil 1^1.2 ((fe,ifc+i)) < c'fc2-f^-3 
and, similarly 

(33) ll(^)^(x, ■)\\w~s.2^^k,k+i)) < ■)\\w^'H(k,k+i)) < c'k'""-'. 

Assume that ll^^ll^^''.2((j^. ^j^.^) = 1- Tiie duality (4), estimates (32) and (33) show 
that 

^max^ ||(^)"((i?- (/f(2iJ- 1)A;2^-V),^) 11^00(0,1) < c'fc^^-^ 
This especially implies that 

II (^R - {H{2H - l)fc2^- V) : Wo'\k, k + 1) ^ W^'\k, A; + 1)|| < C2k'^"-\ 

Let us denote by / : VF^'^((0, 1)) W~^''^{{0, 1)) the natural imbedding. According 
to Lemma 1 (iv) we have ||/||i^5 < oo. We finally obtain 

II {R - {H{2H - I)fc2^-V) : W^^'^ik, k + l)^ W-''\k, k + 1)\\hs 

< \\I\\hs\\{R- {H{2H -l)k^"-^V) ■.W^'\k,k + l)^W^'\k,k + l)\\ 

< C3/c2^-3. 

This establishes (31) and completes the proof of the theorem. □ 

Remark 5. A closer inspection of the above proof reveals that the constant rn in 
Theorem 4 satisfies vh = H\2H — 1| ||x(o,i) ||^h-i/2 2(o i) ' Especially, rn tends to 
zero as if ^ 1/2. Moreover, one also checks that it is possible to choose Sh = 
min(l,2 - 2H). 
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After Theorem 4 it is natural to ask whether similar phenomena take place if only 
one of the intervals in consideration tends to a point. The answer is positive again. 
Heuristically one might expect that the speed of convergence is only half of what it 
was before, and this actually turns out to be true. 

Theorem 6. Let t > 0. Then there are constants r'jj >Q {with r'fj > for H ^ 1/2) 
and > Q such that as e ^ one has 

\l-H I r^r^2-H\ 



cos{<{E^_^^o),Et,e)) = r^(£A)^-^ +0(£'-^) and 

1 
2 



Proof. As in the proof of Theorem 4 we apply scaling, Corollaries 1 and 2, and 
Lemma l(iv) to the effect that it is equivalent to verify in the limit k —>■ oo that we 
have 

(34) =r^A;^-i + 0(A;^-2) and \\S\\hs = r'ak""-^ + Oik""-^). 

Here S = GR, where G stands for the natural isomorphism G : W^^''^{k,k + !)—»• 
WQ''^{{k,k + 1)) provided by (5), s := 1/2 - H, and 

R : Wo''^((-oo, 0)) ^ W-'''^{k, k + 1) 

is the integral operator 

R(f){x) := H{2H -1) .^_2h dy, lor x e {k,k + 1). 

J ~oo ~ y\ 

This time we consider the auxiliary operator V : VFq'^((— oo, 0)) —>■ VF~^'^(fe, k+1), 
where 

V^{x) := H{2H - 1) r '^^fj Off dy, for xe{k,k + 1). 

Thus V is one-dimensional since its image contains only constant functions. 
According to Lemma 2 it holds that 

II !«- I lliy- = -2((-OC,0)) ~ ^'^ 

Hence, by one-dimensionality and the duality (4) we infer that 

\\V : Wo''^((-oo,0)) ^ W-''\k,k + 1)11 = c'k^-\ 

By using again the decomposition S = GV + G{R — V) we deduce, as in the proof 

of Theorem 4, that the one-dimensionality of V and the Hilbert-Schmidt property 
of the natural imbedding W'^^^{{k, k + I)) ^ W-'^'^dk, k + 1)) (where the Hilbert- 
Schmidt norm is independent of k) enable us to deduce (34) as soon as we establish 
that 

(35) \\{R-V) ■.Wq'^{-oo,0) ^W^^^ik,k + 1)\\ <C2k"-^. 

Observe that V — R has the integral kernel u{x,y) := 2{2H — l){\x — — 
\k — Clearly (35) follows from duality and the estimate 

(36) sup ll(^)"u(x,-)||^^-.,2((_^,o))=0(^''~') for a G {0,1}. 



xeik,k+l) 
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In turn, for a = 1 this estimate is a direct consequence of Lemma 3. In order to 
verify it for a = 0, we fix x G (fc, A; + 1) and apply the same Lemma as follows: 



I "''llvK-«.2((-oo,0)) 



J k 

< C II 1^ - ^llvK-s.2((_oo,0))'^* - '^'^^ ^• 

Jk 



In the second inequality above we made use of the Minkowski inequality for Banach 
space norms. □ 

The remaining cases are simpler to handle and they are collected in the following 
theorem. 

Theorem 7. (i) Let H ^ ^. Then /(i?(_e o) '■ E{o,s)) = for any e > 0. 

(ii) Let H Then I{E(^_^^_^) : £^(£,oo)) = for any e > 0. 

(iii) <(-E(-oo,o)>^(o,oo)) > 0- 

(iv) Let t\ <t <t2 he arbitrary. Then for small enough e > it holds that 
(37) -'"(-E'(-oo,ti)u(t2,oo) : Et^s) < ce^~^- 

Proof, (i) Assume the contrary, that is, /(£'(_£ q) ■ ^(o.e)) < °° for some e > 0. 
Since FBM possesses local independence, its infinitesimal space is trivial, that is, 
nj^i -£'(o,±£/n) = {0} (otherwise the Gaussian space would have uncountable dimen- 
sion; see Proposition 5 of [8]). By Theorem 1 of [16], this implies the corresponding 
relation for a-algebras, i.e. H^i c''(-^(o,±£/n)) = up to sets of measure or 1. 

Theorem 1 (iv) then yields that lim„^oo -^(-£'(-£/ri,o) • ^(o.e/n)) — ^{{^S} '■ {^i^}) = 
0. On the other hand, we have I{E(^_^ q^ : -^(ce)) > when i7 / |. Now, however, 
the self-similarity of FBM implies that L{Ei^_^i^ q^ : £'(o^£/n)) does not depend on n, 
and we get a contradiction. 

(ii) By self-similarity. Theorem 1 (iii) and the previous claim, we have 

= -^(£^(-00,0) : -2;(o,oo)) 

> ^(£^(-£,0) : -E^{o,e)) = oo. 

(iii) This is an immediate consequence of Lemma l(i) and Lemma 3, since together 
they imply that for a dense set of elements Xi G -E^(-oo,0) ^2 € -E(o,oo) have 
that 

max(||Xi||,||X2||)<c||Xi-X2||. 

(iv) Write Ai = -^(-octi); ^2 = -E'(t2,-oo); ^ = Et,e- Since the angle between 
the subspaces Ai and A2 is positive, we see that A := span(Ai |J A2) is naturally iso- 
morphic (not necessarily isometric) to the direct sum {AiQ)A2)i2. In this isomorphism 
the operator Pg ■ A B conjugates to the operator [Pb : Ai B, Pb : A2 ^ B], 
whose Hilbert-Schmidt norm is bounded by cs^~^ by Theorem 6. This proves the 
claim. □ 

5. Generalizations and open questions 

The most natural generalization of FBM to R" is the Levy FBM, which is defined 
as the Gaussian process indexed by the parameter u G R"^ and having the 
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covariance structure 

Here H G (0, 1). As in the one-dimensional case this process has a version that has 
Holder continuous realizations. We refer to [6, Chapter 18] for the existence and basic 
properties of n-dimensional Levy FBM. We will sketch the proof of an n-dimensional 
version of Theorem 4. For that end we first present an auxiliary result. 

Lemma 5. Let n > 2 and s G (— n/2 — 1,— n/2). Then there is a constant c > 
such that for every ^ G C^(S(0, 1)) with J^^^dx = and f e C"+i(5(0,l)) it 
holds that 

•^^i^'^> l<\a\<n+l 

Proof. Observe that in the left hand side we may replace / by / — m, where m is 
the average of / over the ball B{Q,1). Hence we may assume that 11/11/^30( ^(0 i) ) 
is dominated by ||-P/|| j^oo( g(o i-^ y It follows that it is enougli to prove the stated 
estimate where one sums over all \a\ < n + 1 in the right hand side. But it is 
easy to extend / to an element / G W^'^^''^{'R^) with norm less than constant 
times Yl\a\<n+i II-^"/IIl°°( b(o 1) )- '^^^ claim follows now by duality since formally 
W+^'^iBiO, 1)) C W-''^{B{0, 1)) = W^'^{B{0, 1))'. □ 

Theorem 8. Let {Xsj^gRn be an n-dimensional Levy FBM with Hurst parameter 
H G (0, 1). For any ball B C R" let Eb be the L'^-space generated by the differences 
{Xg^ —Xs2\si,S2 G B}. Then, if si ^ S2 the subspaces EB[si,e) o-f'-d EB[s2,e) 
asymptotically independent as e ^ 0. Moreover, there are positive constants ci, C2 > 
such that 

Proof. The proof is analogous to the proof of Theorem 4. First of all, the lower 
bound is an immediate consequence of the one-dimensional case since the restriction 
of the process to a line through the points si, S2 is a one-dimensional FBM. In order 

to deduce the upper bound we observe that according to Lemma 4 and an easy 
analogue of Lemma 3 the cosine of the angle between the spaces is given by the 
quantity 



1 



A:=— -sup / / I"" — (l){u)ip{s) duds, 

where the supremum is taken over all functions G C^(-B(0, l))nWo -^'^(^(0, 1)) 
and lb G CQ°{B{kei,l)) n Wq "^^ '^'^(i?(fcei, 1)), with unit norm and zero mean. 
Here k = \si — S2|/e > 0. Observe that we used the obvious scaling and rotation 
invariance of the Levy FBM. By a twofold application of Lemma 5 it follows that 

A< sup \D:D^{\u-sn\^e^-'. 

"e-^(0'l)'^S-^('=^i'l)l<|a|<n+l,l<|/3|<n+l 

□ 

Our results raise several interesting open problems related to local independence 
of stochastic processes. We expect that the methods of the present paper are pretty 
much restricted to dealing with the FBM, although they may help in obtaining 
insights and conjectures regarding the following questions. 
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Q.l Let X = {XijigR, be a Gaussian process with continuous paths and station- 
ary increments. Find necessary and sufficient conditions for the local independence 
property, e.g. in terms of the spectral measure of X, or in terms of the variance 
function v{t) =^X^. 

With regards to Question 1, we can note a couple of obvious obstacles for local 
independence. First, if the process is L^-di£Ferentiable, the value of the derivative 
process belongs to the infinitesimal sigma-algebra around a point (see [16]), and apart 
from trivial cases this will destroy local independence. Second, periodic processes, 
like the periodic Brownian bridge defined by the variance function 

v{t) = = {t mod 1)(1 - {t mod 1)), 

clearly do not satisfy local independence for all times. Periodic components are 
refiected as atoms of the spectral measure. But are non-smoothness and continuity 
of spectrum already sufficient for local independence? 
One can also ask for a local characterization: 

Q.2 Let {Xt) again be a Gaussian process with stationary increments. Give 
conditions on the variance function v{t) in a neighbourhood of the origin and in a 
neighbourhood of the point \ti — t2\ that would guarantee local independence with 
respect to points ti^t^- 

Q.3 Superposing Brownian bridges with different periods, one can probably 
build examples of non-smooth processes where local independence breaks over any 
rational distance. But is it possible to construct a continuous but non-differentiable 
Gaussian process with stationary increments that does not possess local indepen- 
dence over any distance? 

Q.4 So far we have only focused on Gaussian processes. Our information-based 
definition of local independence is, however, meaningful for any kind of stochastic 
process. It is then interesting to ask about the local independence of various depen- 
dent processes. For example, do fractional Levy processes have this property? 
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